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Abstract 

New generalization of ^-scaling in inclusive particle production is proposed. The 



scaling variable z is fractal measure which depends on kinematical characteristics of 
the underlying sub-process expressed in terms of the momentum fractions X\ and £2 
of the incoming protons. In the generalized approach, the X\ and are functions of 
the momentum fractions y a and of the scattered and recoil constituents carried out 
by the inclusive particle and recoil object, respectively. The scaling function ip(z) for 
charged and identified hadrons produced in proton-proton collisions is constructed. 
The fractal dimensions and heat capacity of the produced medium entering definition 
of the z are established to obtain energy, angular and multiplicity independence of 
the ip(z). The scheme allows unique description of data on inclusive cross sections 
of charged particles, pions, kaons, antiprotons and lambdas at high energies. The 
obtained results are of interest to use ^-scaling as a tool for searching for new 
physics phenomena of particle production in high transverse momentum and high 
multiplicity region at proton-proton colliders RHIC and LHC. 
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1 Introduction 



Production of particles with high transverse momenta from collision of hadrons and 
nuclei at sufficiently high energies has relevance to constituent interactions at small scales. 
This regime is of interest to search for new physical phenomena in elementary processes 
such as quark compositeness [I], extra dimensions [2], black holes fractal space-time 
[3] etc. Other aspects of high energy interactions are connected with small momenta 
of secondary particles and high multiplicities. This regime has relevance to collective 
phenomena of particle production. Search for new physics in both regions is one of 
the main goals of investigations at Relativistic Heavy Ion Collider (RHIC) at BNL and 
Large Hadron Collider (LHC) at CERN. Experimental data of particle production can 
give constraints for different theoretical models. Processes with high transverse momenta 
of produced particles are most suitable for precise test of perturbative QCD. The soft 
regime is of interest for verification of non-perturbative QCD and investigation of phase 
transitions in non-Abelian theories. 

Nucleus-nucleus interactions are very complicated. In order to understand their nature 
one often exploits phenomenology and comparison with more simple proton-nucleus and 
proton-proton collisions. Many approaches to description of particle production are used 
to search for regularities reflecting general principles in these systems at high energies 
One of the most basic principles is the self-similarity of hadron production valid 
both in soft and hard physics. Other general principles are locality and fractality which 
can be applied to hard processes at small scales. The locality of hadronic interactions 
follows from numerous experimental and theoretical studies. These investigations have 
shown that interactions of hadrons and nuclei can be described in terms of the interactions 
of their constituents. Fractality in hard processes is specific feature connected with sub- 
structure of the constituents. This includes the self-similarity over wide scale range. 

Fractality of soft processes concerning the multi-particle production was investigated 
by many authors. Fractality in inclusive reactions with high-p^ particles was considered 
for the first time in the framework of z-scaling^S]. The approach is based on principles 
of locality, self-similarity and fractality. It takes into account fractal structure of the 
colliding objects, interaction of their constituents and particle formation. The scaling 
function ip(z) and the variable z is constructed via the experimentally measured inclusive 
cross section Ed 3 a/dp 3 and the multiplicity density dN/drj. In the original version, the 
construction was based on the assumption that gross features of the inclusive particle 
distribution for the inclusive reaction 



at high energies can be described in terms of the corresponding exclusive sub-process 



The Mi and M 2 are masses of the colliding hadrons (or nuclei) and m\ is the mass of 
the inclusive particle. The mass parameter m 2 is introduced in connection with internal 
conservation laws (for isospin, baryon number, strangeness...). The X\ and x 2 are fractions 
of the incoming four-momenta Pi and P 2 of the colliding objects. The scaling variable 
z was constructed as fractal measure in terms of the fractions x\ and x 2 . It depends 
on the nucleon anomalous dimension S and the average multiplicity density dN/dr]\o of 



Mi + M 2 -> mi + X 




(xiMi) + (x 2 M 2 ) -> mi + (xiMi + x 2 M 2 + m 2 ). 



(2) 
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charged particles produced in the central region of the interaction (77 = 0). The scaling 
for non-biased collisions was established for single constant value of 5 as independence of 
the scaling function ip{z) on the collision energy y/s and angle 9 of the inclusive particle. 

The concept of ^-scaling was generalized for various multiplicities of produced particles 
[TKj . Relation of the ^-scaling to entropy S and heat capacity c of the colliding system was 
established. The generalization was connected with introduction of momentum fraction 
y in the final state written in the symbolic way 

(xiMi) + (x 2 M 2 ) -> (m x /y) + (x 1 M 1 + x 2 M 2 + m 2 /y). (3) 

It was shown that generalized scaling represents regularity in both soft and hard regime in 
proton-(anti)proton collisions over a wide range of initial energies and multiplicities of the 
produced particles. However, the generalization of the scaling for various multiplicities 
was at the expense of the angular independence of the scaling function observed for y — 1. 

In this paper we show that independence of the scaling function ip(z) on the collision 
energy y/s, multiplicity density dN/drj and angle 9 can be restored simultaneously, if two 
fractions y a and y& for the scattered constituent and its recoil are introduced, respectively. 

2 New generalization of ^-scaling 

We consider collision of extended objets (hadrons and nuclei) at sufficiently high en- 
ergies as an ensemble of individual interactions of their constituents. The constituents 
are partons in the parton model or quarks and gluons in the theory of QCD. Single in- 
teraction of the constituents is illustrated in Fig.l. Structures of the colliding objects are 
characterized by the anomalous dimensions 5\ and 5 2 in the space of momentum frac- 
tions. Interacting constituents carry the fractions x± and x 2 of the incoming objects. 
The sub-process is considered as a binary collision with production of the scattered and 
recoil constituents, respectively. The inclusive particle carries the momentum fraction 
y a of the scattered constituent which fragmentation is characterized by the anomalous 
dimension e a . Fragmentation of the recoil constituent is described by the corresponding 
fractal dimension e& and momentum fraction y^. 

Multiple interactions are considered to be similar. The property is manifestation of 
the self-similarity of hadronic interaction at the constituent level. The interactions are 
governed by local energy-momentum conservation law. The self-similarity at small scales 
reflects fractality of the interaction objects and their constituents characterized by the 
corresponding fractal dimensions. The fractality concerns parton content of the composite 
structures involved. 

2.1 Locality, self-similarity and fractality 

The idea of 2-scaling is based on the assumption [TU] that gross features of inclusive 
particle distribution of the reaction can be described at high energies in terms of the 
kinematic characteristics of the corresponding constituent sub-processes. We consider the 
sub-process as binary collision 

(xiMO + (x 2 M 2 ) -> (mi/3/ ) + (xiMi + x 2 M 2 + m 2 /y b ) (4) 

of the constituents (x\M\) and (x 2 M 2 ) resulting in the scattered (m\/y a ) and recoil 
(xxMi + x 2 M 2 + m 2 /yb) constituents in the final state. The inclusive particle with the 
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mass mi and the 4-momentum p carries out the fraction y a of the 4-momentum of the 
scattered constituent. Its counterpart (1712), moving in the away side direction, carries 
out the 4-momentum fraction y b of the produced recoil. The binary sub-process satisfies 
the energy-momentum conservation written in the form 

{x 1 P 1 + x 2 P 2 - p/y a f = (xiMx + x 2 M 2 + m 2 /y b f. (5) 

The equation is expression of the locality of hadron interaction at constituent level. It 
represents kinematical constraint on the fractions Xi, x 2 , y a and y b . 

The self-similarity of hadron interactions reflects property that hadron constituents 
and their interactions are similar. This is connected with dropping of certain dimen- 
sional quantities out of description of physical phenomena. The self-similar solutions are 
constructed in terms of self-similarity parameters. We search for the solution 

= TT—^f- (6) 

rK ' Na in dz v ; 

depending on single self-similarity variable z. Here <jj n is the inelastic cross section of the 
reaction and N is the average particle multiplicity. The variable z is specific dimen- 
sionless combination of quantities which characterize particle production in high energy 
inclusive reactions. It depends on momenta and masses of the colliding and inclusive 
particles, structural parameters of the interacting objects and dynamical characteristics 
of the produced system. We define the self-similarity variable z in the form 

z = zqVL- 1 (7) 

where 

n(x 1 ,x 2 ,y a ,y b ) = {l-x^{l-x 2 ) 5 \l-y a ) e *{l-y b y\ (8) 

The variable z has character of a fractal measure. For a given reaction its finite 
part z is proportional to the transverse kinetic energy of the constituent sub-process 
consumed on the production of the inclusive particle (mi) and its counterpart (m 2 ). 
The divergent factor Q^ 1 describes resolution at which the sub-process can be singled 
out of this reaction. The Q(xi,x 2 ,y a ,yb) is relative number of parton configurations 
containing the incoming constituents which carry the fractions X\ and x 2 of the momenta 
Pi and P 2 and the outgoing constituents which fractions y a and y b are carried out by the 
inclusive particle (mi) and its counterpart (m 2 ), respectively. The 5\ (<5 2 ) and e a (e b ) are 
anomalous fractal dimensions of the incoming and fragmenting objects, respectively. For 
inelastic proton-proton collisions we have S\ = S 2 = 5. We also assume that fragmentation 
of the scattered and recoil constituents is governed by the same anomalous dimensions 
e a = e& = e. 

Common property of fractal measures is their divergence with the increasing resolution 

z(Q) ^00 if n- 1 -> 00. (9) 

For the infinite resolution, all momentum fractions become unity {x\ = x 2 = y a = y b = 1) 
and Q = 0. The kinematical limit corresponds to the fractal limit z = 00. 

2.2 Principle of minimal resolution 

The momentum fractions xi, x 2 , y a and y b are determined form principle of minimal 
resolution of the fractal measure z. The principle states that resolution Q' 1 should be 
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minimal with respect to all binary sub-processes (J1J in which the inclusive particle ra\ 
with the momentum p can be produced. This singles out the underlying interaction of 
the constituents. The momentum fractions Xi, x 2 , y a and y b are found from minimization 

of Q' 1 (x 1 ,x 2 ,y a ,y b ), 

dn(xi,yj) dn(xj,yj) dn(xj,yj) dn(xj,yj) 

0x1 ax 2 ay a ay b 

taking into account the energy-momentum conservation (J5J). The momentum fractions X\ 
and x 2 can be decomposed as follows 

a?i = Ai + xi(a), ^2 = A 2 + x 2 (a). (11) 

The parameter a = 82/81 is ratio of the anomalous fractal dimensions of the colliding 
objects. Using the decomposition, the expression (jlj) can be rewritten to the symbolic 
form 

xi + x 2 -> (Ai + A 2 ) + (xi +X2). (12) 

This relation means that A-parts of the interacting constituents contribute to the pro- 
duction of the inclusive particle, while the x _ P ar ts are responsible for the creation of its 
recoil. The X's are functions of y a and y b 

Ai = K X /y a + vi/yb, A 2 = K 2 /y a + u 2 /y b , A = vo/yl - vo/yl, ( i3 ) 

where 

(P2P) = (PiP) (u) 

1 (P 1 P 2 )-MiM 2 ' 2 (PiP 2 )-MiM 2 ' 1 J 

_ M 2 m 2 _ M l m 2 

Vl ~ (P 1 P 2 ) - M X M 2 ' " 2 ~ {P X P 2 ) - M X M 2 ' 1 j 

0.5m? _ 0.5m2 

^ = {P^)-MM U ° = (PiP 2 )-M 1 M 2 - ( } 

The x' s are expressed via X's as follows 



Xi = + -ui, X2 = yi4 + ul + ^2, (17) 

where \ \ I A 

/i 2 = (AiA 2 + A )a- — , /i 2 = (AiA 2 + A )a _1 - r- 2 -, (18) 

1 — A 2 1 — Ai 

and cjj = //jC7 (z = 1, 2). 
The quantity 

U = (19) 
2V« 

has physical meaning of longitudinal component of the " structural 4- velocity" \Tj\ . It is 
function of a and the kinematical factor 



Wi^W^r (20) 
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(0 < £ < 1). The £ characterizes kinematical scale of the underlying constituent interac- 
tion. 

Solution of the system (fTUj) with the condition (jSJ) can be obtained by searching for 
the unbounded maximum of the function 

F(y a ,yb) = tt(x 1 (y a ,y b ),x 2 (y a ,yb),ya,yb) (21) 

of two independent variables y a and y b . Here Xi(y a ,y b ) are given explicitly by the ex- 
pressions (jllj) . There exists single maximum of the function F(y a ,yb) in the allowed 
kinematical region and we determined it numerically. 

2.3 Scaling variable z 

Search for an adequate, physically meaningful but still sufficiently simple form of the 
self-similarity parameter z plays a crucial role in our approach. We define the scaling 
variable z in the form 

1/2 

z= — fT 1 . (22) 

Here m is a mass constant which we fix at the nucleon mass. The transverse kinetic energy 
of the constituent sub-process consumed on the production of the inclusive particle (mi) 
and its counterpart (777.2) is determined by the formula 

sT = T a + T b , (23) 

where 

T a = y a (s\ /2 - M 1 \ 1 - M 2 A 2 ) - m x , (24) 
T b = y b (s]/ 2 - M lX i - M 2X2 ) - m 2 . (25) 
For more details see the Appendix. The terms 



sT = y/(X 1 P 1 + A 2 P 2 ) 2 , s][ 2 = ^{xiPi + X2P2) 2 (26) 

represent energy for production of the scattered constituent and its recoil, respectively. 
The boundaries of the range of the variable z are and 00. They are accessible at any 
collision energy. 

The dN/dr]\o is average multiplicity density of charged particles produced in the central 
region of the reaction ([T]) at pseudorapidity r\ = 0. It depends on state of the produced 
medium in the colliding system. The parameter c characterizes properties of this medium. 
The quantity 

W = (dN/dr]\ ) c ■ VI (27) 

is proportional to all parton and hadron configurations of the colliding system which can 
contribute to production of the inclusive particle with the momentum p. The scaling 
variable (j22j) is proportional to the ratio 

1/2 

- ~ w (28) 

1 /2 

of the transverse kinetic energy s_[ and maximal number of the configurations W. 
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2.4 Scaling variable z and entropy S 

According to statistical physics, entropy of a system is given by number of all statistical 
states W of the system as follows 

S = In W. (29) 
In thermodynamics, entropy for ideal gas is determined by the formula 

S = c v In T + R In V + const. (30) 

The cy is heat capacity and R is universal constant. The temperature T and the volume 
V characterize state of the system. Using (|27jl and (J2HJ), we can write 

S = c\n[dN/d V \ ]+\n[(l-x 1 ) s ^l-x 2 )^(l-y a ) ea (l-y b y b ]. (31) 

Exploiting analogy between Eqs. (fHUj) and (JHH), we interpret the parameter c as "heat ca- 
pacity" of the produced medium. The multiplicity density dN/dn\o has physical meaning 
of "temperature" of the colliding system. The second term in Eq. (|31|) depends on volume 
in space of the momentum fractions {xi, x 2 , y a , Ub}- This analogy emphasizes once more 
interpretation of the parameters 5±, 5 2 , e a and e& as fractal dimensions. As seen from Eq. 
(JHIJ), entropy of the colliding system increases with the multiplicity density and decreases 
with increasing resolution Q^ 1 . 

Let us note that entropy (|29|) of a system is determined up to an arbitrary constant 
lnM^o- Dimensional units entering definition of the entropy can be included within this 
constant. Therefore, it allows us to make analogy between dimensionless dN/dn\ and the 
temperature. This degree of freedom is connected with the transformation 

z -> W ■ z, V -> Wo 1 ■ ip. (32) 

In such a way the scaling variable and the scaling function are determined up to an 
arbitrary multiplicative constant. 



2.5 Scaling function ip{ z ) 

The scaling function ip(z) is expressed in terms of the experimentally measured inclu- 
sive invariant cross section Ed 3 a/dp 3 , multiplicity density dN/dn and the total inelastic 
cross section cr,j n . Exploiting the definition © one can obtain the expression 

its i 1 ^d 3 a 

^ = -(dN/d V )a m J ~ E W (33) 

Here s is the centre-of-mass collision energy squared and 

= drj_dz_ _ dr]_dz_ 
8ki 8k 2 8k 2 8ki 

is the corresponding Jacobian. Angular properties the ip(z) showed that the scaling is 
valid rather in pseudorapidity than in rapidity (see Appendix). The function ip(z) is 
normalized as follows 

ip{z)dz = 1. (35) 
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The relation allows us to interpret the ip{z) as a probability density to produce inclusive 
particle with the corresponding value of the variable z. 

3 Properties of the scaling function ip(z) 

Let us investigate properties of z presentation of experimental data obtained in proton- 
proton collisions at high energies. 

3.1 Energy independence of ip(z) 

We analyze experimental data [111-123 on inclusive hadron (h^, tt~, K~ and p) produc- 
tion in minimum-biased proton-proton collisions. The data on inclusive cross sections were 
measured in central rapidity region at FNAL, ISR and RHIC energies ^/s = 19 — 200 GeV. 

The energy dependence of the charged hadron spectra on the transverse momentum 
is shown in Fig. 2(a). The distributions cover the range up to px — 10 GeV/c. The 
cross sections change within the range of 12 orders of magnitudes. Strong dependence of 
the spectra on the collision energy y/s increases with transverse momentum. Figure 2(b) 
shows z-presentation of the same data. The scaling variable z depends on the average mul- 
tiplicity density of charged particles produced in the central pseudorapidity region of the 
collision. We have used experimentally measured values of dN/drjl^o for minimum- 
biased collisions in the analysis of energy and angular properties of ip{z). Independence 
of the scaling function ip{z) on collision energy y/s is found for the constant values of the 
parameters c = 0.25,5 = 0.5 and e = 0.2. The form of ip(z) manifests two regimes of 
particle production. The hard regime is characterized by the power law ip{z) ~ z~@ at 
large z. Soft processes correspond to behavior of the ip{z) at small z. Slope of the scaling 
curve decreases with z in this region. 

Invariant cross sections for 7r~-meson production as function of the collision energy and 
transverse momentum are plotted in Fig. 3(a). The spectra were measured over a wide 
transverse momentum range px — 0.1 — 10 GeV/c. The cross sections change from 10 2 to 
10~ 10 mb/GeV 2 . Strong dependence of the pion spectra on y/s as for charged hadrons is 
observed. The ^-presentation of the same data is shown in Fig. 3(b). For pions (as well 
as for all other types of the particles - kaons, antiprotons, ... ), the dependence of z on 
the charged particle multiplicity density dN/dr]\ v=Q have been used in the formula (J2*2*|) . 
Let us stress that, unlike this, the scaling function (|3*3*|) is normalized to the multiplicity 
density of pions. Independence of the scaling function for pions on y/s was obtained at 
c = 0.25, 5 = 0.5 and e = 0.2 as for charged hadrons. Shape of the ip(z) is similar in both 
cases, as well. 

Transverse momentum spectra for i^~-meson production are shown in Fig. 4(a). The 
cross sections were measured in the range px = 0.1 — 8 GeV/c. Data for if^-mesons 
obtained by the STAR Collaboration at RHIC are presented in the Fig. 4(a) as well. 
The i^-meson spectra demonstrate strong dependence on the collision energy y/s. The 
corresponding scaling function ip(z) is depicted in Fig. 4(b). Independence of the ip(z) on 
y/s is restored at c = 0.25, 5 = 0.5 and e = 0.3. Similar features of px and z presentations 
of experimental data on antiproton production are presented in Figs. 5(a) and 5(b). 
The energy independence of if)(z) for antiprotons is established at c = 0.25,5 = 0.5 and 
e = 0.35. 
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As a result we can conclude that energy independence of the scaling function ip{z) 
is valid for different types of hadrons in a wide range of centre-of-mass energy y/s and 
transverse momentum p?. 

3.2 Angular independence of ij){z) 

We analyze experimental data ^01 GZj on angular dependence of negative hadrons 
(pions, kaons and antiprotons) measured at ISR energies. The data were measured both 
in the central and fragmentation regions. Results of our analysis is demonstrated at the 
energy yfs = 53 GeV. 

Invariant cross sections for 7r~-meson production as function of the centre-of-mass 
angle 9 and transverse momentum are shown in Fig. 6(a). The angles cover the range 
9 = 3°— 90°. The central and fragmentation regions are distinguished by different behavior 
of differential cross sections. The z-presentation of the same data is demonstrated in Fig. 
6(b). The charged hadron multiplicity density dN/dr]\ v=0 represents angular independent 
factor in the definition of the variable z. Contrary to this, the scaling function (}3*3*j) is 
normalized to the multiplicity density dN/dr] of pions depending on the angle 9. Angular 
and energy independence of the scaling function for pions was obtained at the same values 
of c = 0.25, 5 = 0.5 and e = 0.2. The function if)(z) for small 9 is sensitive to the value of 
m 2 . This parameter is determined from the corresponding exclusive reaction 

p + p -y tt~ +p + A ++ . (36) 

The reaction is limiting case of the sub-process (j3J) for x\ — X2 — y a — Vb — 1- According 
to Eqs. (0) and (JHHJ), we obtain m 2 = m(A ++ ) — m(p) = 0.3 GeV. This value was used 
in our analysis for inclusive 7r~-meson production. 

Transverse momentum spectra for i^-mesons and antiprotons produced in pp collisions 
at different angles are shown Fig. 7(a) and 8(a). In addition to ISR data at y/s = 53 GeV, 
the data from RHIC at y^s = 200 GeV are shown as well. The angular dependence of 
the spectra demonstrate strong difference between central and fragmentation regions. 
The corresponding function ip{z) for kaons and antiprotons is plotted in Figs. 7(b) and 
8(b), respectively. For both particles, the charged hadron multiplicity density dN/drj\ v=0 
represents angular independent factor in the definition of the variable z. The scaling 
function (}3*3*j) is normalized to the angular dependent multiplicity density dN / drj of kaons 
and antiprotons, respectively. 

Independence of the ip{z) on the angle 9 is obtained at c = 0.25, 5 = 0.5 and e = 0.3 for 
kaons and c = 0.25, 5 = 0.5 and e = 0.35 for antiprotons. Values of these parameters allow 
us to obtain simultaneously angular and energy independence of the scaling function. Like 
in the case of 7r~-mesons, the function ip{z) for kaons and antiprotons is sensitive to the 
value of m 2 at small angles 9. The corresponding exclusive reactions 

p + p^K~+p + p + K + , (37) 

p + p^p + p + p + p (38) 

were used to determine the parameter m 2 . Exploiting Eq. © for x\ = x 2 = y a = Vb = 1 
and Eqs. (|3*7)l and (|3*8|). we obtain m 2 = m(K + ) = 0.5 GeV and m 2 = m{p) = 0.94 GeV 
for the inclusive production of i^~-mesons and antiprotons, respectively. These values 
were used in our analysis. 
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3.3 Multiplicity independence of ij)(z) 

The STAR Collaboration obtained the new data [28J on multiplicity dependence of 
the inclusive spectra of charged hadrons produced in pp collisions in the central rapidity 
range \rj\ < 0.5 at the energy y/s = 200 GeV. The transverse momentum distributions 
were measured up to 9.5 GeV/c using different multiplicity selection criteria. Figure 9(a) 
demonstrates strong dependence of the spectra on multiplicity density at dN/dr] = 2.5, 6.0 
and 8.0. The same data are presented in Figure 9(b) in the scaling form. The scaling 
function ip(z) changes over 6 orders of magnitude in the range of z = 0.2 — 10. The 
independence of ip(z) on multiplicity density dN/dr] is obtained. The result gives strong 
restriction on the parameter c. It was found to be c = 0.25. 

The STAR Collaboration measured the multiplicity dependence of the i^-meson and 
A s -baryon spectra [2H] in the central rapidity range \rj\ < 0.5 at the energy y/s = 200 GeV. 
The spectra are presented in Figs. 10(a) and 11(a). The multiplicity density was varied 
in the range dN/dr] = 1.3 — 9.0. The transverse momentum distributions were measured 
up to 4.5 GeV/c. The corresponding function ip(z) is shown in Figs. 10(b) and 11(b), 
respectively. The scaling behavior gives strong restriction on the value of c. Data prefer 
c = 0.25 in both cases. 

Thus we conclude that the available experimental data on the multiplicity dependence 
of spectra of charged hadrons, i^-mesons and A s -baryons produced in pp collisions at 
RHIC confirm generalized ^-scaling for the same value of the parameter c. 

4 Conclusions 

Generalized z-scaling for the inclusive particle production in proton-proton collisions 
was suggested. The scaling variable z is function of multiplicity density dN / dr)\o of charged 
particles in the central region of collision. The variable z depends on the parameters c, 5 
and e. They are interpreted as specific heat of the produced medium, anomalous fractal 
dimension of the proton and fractal dimension of the fragmentation process, respectively. 
Connection between the scaling variable z and entropy S of the interacting system was 
established. 

We have analyzed experimental data on inclusive cross sections of hadrons (h , 7r~, 
K~ , Kg, p and A) measured in proton-proton collisions at FNAL, ISR and RHIC. The 
data cover a wide range of collision energy, transverse momenta and angles of the pro- 
duced particles. Spectra from minimum biased events and events with various multiplicity 
selection criteria have been studied. The energy, angular and multiplicity independence 
of the scaling function was established. It gives strong constrains on the values of the 
parameters c, 5 and e. It was shown that the parameters are constants in the considered 
kinematical region. The parameters c = 0.25 and 5 = 0.5 were found to be the same for 
all types of the considered inclusive hadrons. The value of e increases with mass of the 
produced hadron. 

The variable z has property of a fractal measure connected with parton content of the 
composite structures involved. The fractal dimensions S and e determine fractal properties 
of z in space of the momentum fractions. The scaling function ip{z) manifests two regimes 
of particle production. The hard regime is characterized by the power law ip{z) ~ at 
large z. Soft processes correspond to behavior of the ip{z) at small z. Slope of the scaling 
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curve decreases with z in this region. On basis of the performed analysis we conclude that 
^-scaling in proton-proton collisions is regularity which reflects self-similarity, locality and 
fractality of hadron interaction at constituent level. It concerns structure of the colliding 
objects, interactions of their constituents and fragmentation process. 

We consider that obtained results are of interest for searching and study of new physics 
phenomena in particle production over a wide range of collision energies, high transverse 
momenta and large multiplicities in proton-proton and nucleus-nucleus interactions at the 
RHIC and LHC. 

Acknowledgments. The investigations have been partially supported by the 
IRP AVOZ10480505, by the Grant Agency of the Czech Republic under the contract No. 
202/04/0793 and by the special program of the Ministry of Science and Education of the 
Russian Federation, grant RNP.2. 1.1.5409. 



5 Appendix 



The invariant differential cross section for production of the inclusive particle is nor- 
malized as follows 

d 3 a 2 

E—dr]d p ± = a ind N. (39) 

The (T ine i is the inelastic cross section and N is the average multiplicity. The differential 
cross section of identified hadrons of certain type is normalized to the corresponding 
average multiplicity of this type. The inclusive cross section can be expressed in terms of 
K\ and k.2 in the way 



d 3 a 1 J(PiP 2 ) 2 -M?M$ d 2 a 

E = . (40) 

dp 3 2tt [(PiP 2 ) - M X M 2 ] 2 dn x dK 2 v ' 

In the region of high energies, the formula can be written in the approximate form 

d 3 a 1 d 2 a , . 

dp 6 us dKidK 2 

where s is square of the centre-of-mass energy. We suppose that the inclusive cross section 
is given by solution (jOJ) as function of a single variable z = z(ki, k 2 ). Another independent 
combination of K\ and k 2 is rapidity 

y = \^-. (42) 

2 K\ 

Using the variables z and y, we get the normalization 

d 2 cr f d 2 <j f 

-dK X dK 2 = / ——dydz = o inel \ p(y)ip(z)dydz = a ine iN, (43) 



dnidn 2 J dydz 

where p(y) = dN/dy is the rapidity distribution of particles of the considered type. 
Detailed analysis of experimental data on angular properties of the scaling function showed 
that the factorization 

Si = a *" wmiz) (44) 



ii 



is valid using rather pseudorapidity than rapidity. In the centre-of-mass of the symmetric 
systems (e.g. for NN collisions) we can exploit the relations 

1 1-COS02 1 K 2 

rt = -In — ~ y — -In — , (45) 

' 2 l + cos6 l y 2 kx K ' 

_ (P 2 p) Ep + pcosOx _ (J\p) E p -p cos fl 2 / 4fi x 

Kl ~ (P1P2) - M X M 2 ~ ^ ' K2 ~ (PiPg) - M X M 2 ^ 1 j 

and set at the end X = 2 = 0. Then we get 

<9?7 <9?7 (90i >/i 1 



9«i d9\dK\ 2p (l+cos6) 



(47) 



<9k 2 90 2 9k 2 2p (1 -cos0)' 1 ' 

The derivatives dz/di^i which enter the Jacobian (|34j) were calculated numerically. 

1 /2 _____ 

Finally, we examine the transverse kinetic energy s_[ which enter the definition f!22l) of 
the scaling variable z. The transverse momentum balance in the constituent sub-process 
is guarantied by the identity 

X1X2 = A*iA*2 = AiA 2 + A . (49) 

Neglecting masses Mi\ and MiXi of the interacting constituents in the expressions (J24|) 
and (|25|). we get for the transverse kinetic energies T a and the following relations 

T± ~ y aV /AiA 2 2P 1 P 2 - mi ~ ^(Pi-) 2 + m i - mi ( 5 °) 

and 

Tj ~ y b J X1X22P1P2 - m 2 ~ y b J(\i\ 2 + A )s - m 2 



Vb\j KiK 2 s/y 2 a + m 2 2 /yl - m\/y 2 a - m 2 



~ 2/6 V ((Pi) 2 + m l)/?/a + - m lA/a ~ m 2 = \j G»±) 2 + ™% - 77l 2 . (51) 

In the last equation we have used the transverse momentum balance 

Pl/Va=p h ±/Vb- (52) 

valid in the constituent sub-process. The p a and p b are momenta of the inclusive particle 
(mi) and its counterpart (m 2 ), respectively. 
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Figure 1. Diagram of the constituent sub-process. 
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(b) The corresponding scaling function ip(z). 
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